Operational regimes of lasers based on gain media with a large Raman scattering cross-section

We report on unusual regimes of operation of a laser with a gain medium with a large Raman scattering cross-section, which is often inherent in new types of gain media such as colloidal and epitaxial quantum dots and perovskite materials. These media are characterized by a strong electron–phonon coupling. Using the Fröhlich Hamiltonian to describe the electron–phonon coupling in such media, we analyze the operation of the system above the lasing threshold. We show that below a critical value of the Fröhlich constant, the laser can only operate in the conventional regime: namely, there are coherent cavity photons but no coherent phonons. Above the critical value, a new pump rate threshold appears. Above this threshold, either joint self-oscillations of coherent phonons in the gain medium and photons in a cavity or a chaotic regime are established. We also find a range of the values of the Fröhlich constant, the pump rate, and the resonator eigenfrequency, in which more than one dynamical regime of the system is stable. In this case the laser dynamics is determined by the initial values of the resonator field, the active medium polarization, the population inversion, and phonon amplitude.

Recent advances in chemistry and nanotechnology have significantly expanded the types of active laser media. At present, active media based on colloidal 1,2 and epitaxial quantum dots (QDs) [3][4][5] , perovskites [6][7][8] , and dye molecules [9][10][11] are characterized by a high gain and are widely used for the miniaturization of lasers [12][13][14][15][16] . As a rule, the interaction between electrons and optical phonons is strong in active media with large Raman crosssection [17][18][19] . In this case, it becomes possible to amplify phonons through the action of the electromagnetic field on the electronic subsystem. Thus, in addition to high light amplification, systems with a strong electron-phonon coupling can effectively amplify optical phonons with a short wavelength and create phonon nanolasers 20 . This is facilitated by the fact that the wavelength of an optical phonon at terahertz frequencies is equal to several nanometers, which corresponds to the typical length scale of nanostructures 21 . Since phonons have a relatively high mean free path (of 100-1000 μm) 22 in the terahertz range, they can provide high resolution when used for image processing and nanostructure research.
In this work, we focus on the influence of nuclear vibrations in materials with a strong electron-phonon coupling upon the amplifying properties of conventional lasers. To be specific, we consider amplifying media based on QDs. In such lasers, the pumping is assumed to be incoherent, and lasing is determined by the interaction of the polarization of the QD electronic subsystem with the resonator mode. The coupling of the cavity electric field and the active medium polarization is treated in the dipole approximation and is characterized by the Rabi constant 23,24 . As a rule, the observed gain in conventional lasers is much greater than the gain due to the stimulated Raman scattering (SRS). Therefore, even in lasers using active media with a strong interaction between electrons and phonons, the influence of optical phonons on laser operation is usually neglected.
We describe the interaction between the electrons and optical phonons via the Fröhlich Hamiltonian 18,[25][26][27][28] . Because of the strong coupling of optical phonons with the electric field, the Fröhlich constant of the interaction, g, may be comparable or even larger than the Rabi constant.
We use computer simulation to show that there is a region bounded by the curve g = g cr (D 0 ) in the plane of parameters g, D 0 , where D 0 is the population inversion characterizing the incoherent pump. Inside this region, only usual lasing regimes are possible-these are the trivial regime without the coherent optical field and the lasing regime, in which the mean coherent optical field is not zero. In both regimes, coherent phonons are not excited.
We focus our attention on the regions in which D 0 is greater than the laser threshold D th_opt and the Rabi frequency is much less than optical frequencies. Outside of the region g < g cr (D 0 ) , new operational regimes appear. One of these new regimes is characterized by the excitation of coherent joint self-oscillations of both an optical field in the cavity and optical phonons in the gain medium.
We also find that for QDs with a low transverse relaxation rate, a regime with coherent phonons may be transformed into a chaotic regime with a spectrum containing an infinite number of frequencies incommensurable with optical phonon frequencies.
Finally, we show that in contrast to the conventional lasers, there exists a region of values of g and D 0 values in which two different regimes are stable simultaneously. These regimes may be any pair of the following: the trivial non-lasing regime, the conventional laser regime, the lasing regime accompanied by coherent phonons, and the chaotic regime. Regime which is realized in practice is determined by the initial conditions.

Model.
We study the dynamics of a laser consisting of a single-mode cavity with the frequency ω a and QDs serving as an active medium. In the QD, the electronic degree of freedom is usually an exciton excitation. The exciton interacts with the quantized vibrations of the QD nuclei-optical phonons; we consider one optical phonon mode with the eigenfrequency ω b . Since underlying lasing is a resonant phenomenon, we can neglect the excitation of off-resonant levels and only consider two electronic states of the active medium 29 -the ground state g , which is in the valence band, and the excited state |e� , which is in the conduction band. The transition frequency ω σ between these states corresponds to the excitation of an exciton.
In colloidal QDs and dyes in the optical range, the transverse relaxation rate γ σ is about 0.001-0.01 eV [30][31][32] . This is usually much greater than the decay rate in the resonator γ a . The relationship between the relaxation rates γ a ≪ γ σ takes place in gas lasers 29,33 . We, therefore, can use the gas-laser model to describe the QD-based lasers.
We emphasize that we consider a conventional laser scheme that differs from the Raman laser. As opposed to the case considered here, in the latter, pumping is coherent, and the laser must have at least two cavity modes with the frequency difference �ω equal to the phonon frequency ω b . Because the nonlinear Fröhlich interaction provides the energy flow from the mode with higher frequency to the mode with lower frequency, one does not have to create population inversion in the active medium.
Below we consider a laser whose cavity size L is small enough to ensure inequality c/L ≫ ω b . Since �ω ∼ c/L , this inequality provides �ω ≫ ω b , and therefore, the effect of Raman laser can be neglected. In this case, only one cavity mode is involved in the laser dynamics, and population inversion is necessary. Below we assume that Before proceeding to the study of the effect of the electron-phonon interaction on the operation of a laser, we recall the main points of the theory of lasers 24,34 and introduce the corresponding notations.
The system Hamiltonian is the usual Jaynes-Cummings Hamiltonian with two additional terms describing the energy of intermolecular vibrations of nuclei and the Fröhlich Hamiltonian of the interaction of the electronic and vibrational subsystems 18,25-28 : In Eq. (1), the first term is the Hamiltonian of the electromagnetic field in the cavity, the operators â † and â are the operators of creation and annihilation of a quantum of the electromagnetic field. These operators satisfy the commutation relation â,â † =1 . The second term describes the Hamiltonian of the QD exciton modeled by a two-level system. The operators σ = g �e| and σ † = |e� g are transitions operators from the excited state to the ground state and vice versa, and σ †σ is the operator of the population of the excited electronic state in the QD. These operators satisfy the commutation relation σ † ,σ =D , where the operator D determines the difference between the populations of the excited and ground states. The third term is the interaction of the QD with the cavity mode 24 , � R = −d eg · E(r)/ is the Rabi constant of the interaction, d eg is the matrix element of exciton transition, and E(r) is the electric field "per one photon" of the cavity mode at the QD location r . The fourth term in Hamiltonian describes the optical phonon in the harmonic approximation, b † and b are the operators of creation and annihilation of the phonon; these operators satisfy the commutation relation b ,b † =1 . The last term is the interaction of the electronic and vibrational subsystems 18,[25][26][27][28] . The operator b † +b has the meaning of the operator of the amplitude of the nuclear excitation. The Fröhlich constant can be evaluated employing the expression for the Raman scattering cross-section at the Stokes frequency 18 : The experimentally measured value of σ Cross_ sec tion for CdSe QDs is σ Cross_ sec tion ≈ 4 · 10 −6 Å 2 (see Ref. 19 ).
Using the values of the phonon frequency ω b ≈ 0.025 eV (which corresponds to the optical phonon in CdSe 35 ) and the QD dipole moment d eg ≈ 5 − 10 D , we obtain that the Fröhlich constant is of the order of g ≃ 10 −2 eV which we use in our study.
Thus, Hamiltonian (1) is applicable for systems consisting of active medium with high Raman scattering cross-section placed in the cavity. An example of such an active medium is CdSe, CdS, and PbS quantum dots. Note that, for simplicity, we consider only one optical phonon mode in an active medium.
(2) σ Cross_ sec tion = 4π 3 www.nature.com/scientificreports/ To describe relaxation processes, it is necessary to introduce reservoirs 36,37 which interact with electronic and phonon subsystems of the QD and the cavity mode. We consider reservoirs describing Ohmic losses in the cavity, energy relaxation of the exciton, energy relaxation of optical phonons, and the phonon reservoir responsible for exciton dephasing. Eliminating the reservoir degrees of freedom in the Born-Markov approximation results in the master equation in the Lindblad form 36,37 (for details, see Master equation for the system density matrix in "Methods" section).
Employing Lindblad equation, one can obtain the equation of motion for the expected values of the involved operators. To find the operation regimes, we study the evolutions of the expected values of the following operators: the annihilation operator of the field in the resonator â = a , which is proportional to the average field in the resonator, the operator of the exciton total dipole moment, σ = σ , the operator of the population inversion D = D of the excited level, and the annihilation operator of optical phonon b = b . As a result, we arrive at the following equations of motion, which are the expanded Maxwell-Bloch equations (see Derivation of Eqs. (3)-(6) in "Methods" section): The quantities γ a , γ D , and γ p have the meaning of the rates of respective relaxation and pump processes, γ σ = γ D + γ p + 2γ deph where γ deph is the dephasing rate 37 . In Eq. (6), D 0 = (γ p − γ D )/(γ p + γ D ) . We have neglected the correlations between the operators â and σ and the operators â and D , assuming that âD ≃ â D and â †σ ≃ â † σ . In addition, we have neglected the correlation between operators b and σ assuming that b +b † σ ≃ b +b † σ . By splitting the correlators, we neglect quantum fluctuations; in particular, we do not explicitly take into account spontaneous emission. This is correct if we consider the laser dynamics well below or well above the threshold 18,23,34,36 . Solving Eqs. (3)-(6) does not allow for a detailed description of the spontaneous emission process 38 but describes rather well the induced emission 23,37 . Note that equations (3)-(6) describe an indirect coupling between the cavity electric field of the cavity mode and the optical phonon of the active medium. Indeed, Eq. (3) includes the interaction of a with σ through the term −i� R σ , while Eq. (4) has the term −ig(b + b * )σ that couples σ with b . The direct interaction is absent because we assume that QDs are Raman active.
Below, we describe different operational regimes of a laser by three parameters of the problem. The first parameter is the Fröhlich interaction constant g . At g = 0 , we deal with a conventional laser. The second relevant parameter is the pump rate γ p . This parameter determines the transition to the lasing regime. As commonly accepted, we describe the pump rate with the parameter D 0 = γ p − γ D / γ p + γ D , which is uniquely determined by γ p . The third parameter of our problem is resonator eigenfrequency ω a .

Results
The regime of convention laser without excitation of coherent phonons. The trivial solution that corresponds to the absence of coherent oscillations of a, σ , D , and b has the form This solution is analogous to the trivial solution for a laser, in which the phonon degree of freedom of the active medium molecule is not taken into account. Usually, this solution is stable below the lasing threshold. In the case under consideration, the trivial solution is characterized by a constant expected value of the phonon operator b , which corresponds to a constant displacement of the nuclei.
By increasing the pumping D 0 , we expect, besides the trivial solution, to find a self-oscillating solution. Our computer simulation shows no coherent phonons in this regime, and the coherent part of b , which we denote as b ph , is still equal to zero.
Since we are interested in the stationary behavior of the system, namely, the behavior at large times, , we search for self-oscillating solutions in the trial form: www.nature.com/scientificreports/ with still unknown frequency ω gen and the laser threshold D th_opt . Here a opt and σ opt are time-independent amplitudes of oscillations of a and σ . Thus, according to Eqs. (9) and (13) Substituting Eqs. (11)-(13) into Eqs. (3)-(6) and equating the real and imaginary parts of Eq. (4) to zero, we obtain the following expression for the population inversion D th_opt and the self-oscillation frequency ω gen : where � gen = ω gen − ω a and δ gen = ω gen − ω σ . Equations (16) and (17) have terms proportional to g 2 , resulting in the dependence of the stationary population inversion and the frequency of self-oscillations on the Fröhlich constant g. For the field amplitude in the resonator and the polarization, we obtain the following expressions: The quantity D th_opt has the meaning of the generation threshold for self-oscillation of the cavity field. Such a behavior corresponds to the Hopf bifurcation 39 .
The laser regime with excitation of coherent phonons. In this section, we show that the interaction of an exciton with the optical phonon may lead to the appearance of either a regime with coherent phonons or deterministic chaos.
The resonant case, ω a = ω (g) σ . Let us first consider the resonant case. For this purpose we choose the value of ω a to be equal to the effective exciton frequency ω (g) σ , defined in Eq. (18), for any value of g and D 0 > D th_opt . In this case, we have two parameters g and D 0 that determine operational regimes.
The possible dynamical regimes found with the aid of computer simulations of Eqs. (3)-(6) are shown on plane g, D 0 , in Fig. 1. It turns out that there are two threshold values of pumping D 0 , which separate different regimes. When the pumping D 0 is less than the optical threshold D th_opt determined by Eq. (16), computer simulation shows that only trivial solution Eqs. (7)-(10) exists and is stable (the blue area in Fig. 1). Note that at large pumping, disruption of generation occurs. This is due to the increase in the dissipation rate of the exciton polarization γ σ with an increase in the pump rate γ p . When the pump rate becomes as high as γ p ≫ γ D , γ deph , the dissipation of the polarization becomes so large that laser generation is hindered 33,40 . www.nature.com/scientificreports/ When D 0 > D th_opt , numerical analysis shows of Eqs. (3)-(6) that the trivial solution becomes unstable, and a new self-oscillating solution without coherent phonons arises (see Eq. 20). This solution is stable for any g until D 0 < D th_ph g = ∞ (the yellow region in Fig. 1). In this regime, the expected value of the nuclear displacement due to their coherent vibrations about equilibrium positions is zero, b ph = 0 . Therefore, despite accounting for the electron-phonon interaction, the laser radiation spectrum contains neither Stokes nor anti-Stokes frequencies.
It means that nuclei do not take part in self-oscillation. This is exactly what is usually assumed when neglecting the electron-phonon interaction.
In the region D 0 > D th_ph g = ∞ , our computer simulation of Eqs. (3)- (6) shows that there is a curve D 0 = D th_ph g of the Hopf bifurcations (the solid curve in Fig. 1). This curve bounds the area (colored in orange in Fig. 1) in which the lasing regime without coherent phonons becomes unstable, and as we show below, a new regime with coherent phonons arises. We also show that this is a new Hopf bifurcation leading to the formation of new self-oscillation at the frequency close to ω b .
To study the latter regime, it is convenient to present the trial form for the expected values of the operators â , σ , D , and b in the following form: Here, we split the amplitudes of a(t), σ (t), D(t) and b(t) into the time-independent parts ( a opt , σ opt , D th_opt , and −g D th_opt + 1 /(2ω b − iγ b ) ), corresponding to the conventional laser regime (see Eqs. (11)-(14)), and the time-dependent parts ( a ph (t), σ ph (t), D ph (t) and b ph (t) ), forming new non-harmonic self-oscillations with the period T = 2π/ω b,gen . The lowest frequency ω b,gen is the frequency of the excited coherent phonons ω b,gen ≡ 2π/T . Note that ω b,gen = ω b .
We solve Eqs. (3)-(6) numerically for different values of g and D 0 and consider the presence of time-dependent parts mentioned above.
In order to characterize the transition to the self-oscillation with coherent phonons, we average the quantities a ph (t) , σ ph (t) , |D th (t)| , and |b th (t)| over the period of the self-oscillations: �...� = 1 is the time during which the system reaches stationary oscillations, T is the period of new self-oscillations. Figure 2 shows the dependence of b ph on the pump rate D 0 for fixed g > g cr . This dependence behaves like ∼ D 0 − D th_ph , which shows that the second Hopf bifurcation occurs at the second threshold D th_ph .  www.nature.com/scientificreports/ In the new regime, the spectrum of the self-oscillations is not characterized by a single frequency ω b,gen ; the periodic oscillations are not harmonic. These oscillations can be presented as the Fourier series a ph (t) = n a n exp −iω b,gen nt , where the coefficients a n are time-independent. Figure 3a shows the dependence of the stationary self-oscillations of a ph on time t for eight periods. Note that we use the term "slowlyvarying amplitude" for the complex variable a opt + a ph (t) , which has the absolute value and the phase. This results in the modulation of the total phase as well as the phase of phonon part a ph (t) (see Fig. 3b). The period of the modulation is the same as the period of the modulation of the amplitude of a ph (t) . The reason for this matching is that both oscillations are due to auto-oscillations of the optical phonons and the Fröhlich interaction between phonons and active medium polarization. Figure 4a shows the spectrum of the oscillations of a ph (t) which is determined by the equality Note that this spectrum determines the spectrum of laser radiation 37,41 . The spectrum contains frequencies that are multiples of the frequency of coherent phonons ω b,gen , which is a consequence of the interdependence between a and b . The latter is due to the nonlinearity of the Frohlich interaction between the optical phonon and the exciton. Figure 4b shows the projection of the closed cycle of the self-oscillations on the space |a| , |b| , and D . One can see that the system dynamics is presented by a closed curve which is the manifestation of regime periodicity.
Note that as shown on the plane g, D 0 (Fig. 1), the laser threshold D th_opt and the one for self-oscillations with coherent phonons, D th_ph g , are different. This means that the excitation of coherent phonons is not a direct consequence of the Raman scattering of the coherent cavity mode.
The appearance of b ph = 0 implies the emergence of coherent phonons. The latter does not appear in spontaneous Raman scattering even if the incident field is coherent. Thus, this fact cannot be attributed to the Raman scattering of the laser mode. This is confirmed by our computer experiment. Figure 1 shows that the threshold for a coherent harmonic laser mode D th_opt and the threshold D th_ph for the appearance of coherent phonons are different.  www.nature.com/scientificreports/ Usually, the coherent phonons arise in systems at the coherent anti-Stokes Raman spectroscopy (CARS) or after the onset of Raman lasing. In CARS, the appearance of coherent phonons is associated with the resonant excitation of optical phonons at the frequency difference of two incident coherent waves 18,42,43 . In a conventional laser and in our scheme, the pumping is incoherent, and the only coherent wave is the resonator mode above the laser threshold D th_opt . This is clearly not enough for the CARS mechanism to lead to the excitation of coherent phonons.
The case of detuning ω a − ω (g) σ > 0. As noted above, at g = 0 , instead of ω σ , the frequency of selfoscillations is determined by ω (g) σ , which depends on the values of g , D th_opt , and ω a (see Eq. 18). In this subsection, we consider the operational regimes in the case of ω a = ω (g) σ . For conventional lasers ( g = 0 ), the optimal condition for lasing is matching the resonator frequency ω a with the transition frequency ω σ . When this happens, the laser threshold is minimal. With an increase in the difference |ω a − ω σ | , the lasing threshold D th_opt increases. There is the value of |ω a − ω σ | , at which D th_opt becomes greater than D 0 . Then, the pump rate is not sufficient to start lasing, and self-oscillation does not appear. The trivial solution becomes stable again; consequently, the lasing does not occur.
As our computer simulation shows, similar behavior is observed at g = 0 . In this case, however, the role of ω σ is played by ω (g) σ . At a fixed pump rate D 0 > D th_opt , with an increase in the difference ω a − ω (g) σ , the lasing threshold also grows, and sooner or later, this should lead to the suppression of lasing. By computer simulation, we determine for each g the interval of frequencies ω a , where lasing is observed. The lasing interval for a conventional laser now transforms into the lasing corridor on the plane ω a , g (see Fig. 5).
The results of the computer simulation are shown in Fig. 5a. On the plane of the parameters ω a and g , the dependence ω (g) σ g at fixed D 0 > D th_opt is shown by the blue solid line. The corridor of parameters for which the lasing is observed (the lasing corridor) is bounded by the dotted lines. Outside of this corridor, only the trivial solution is stable, and the system does not lase. Below we deal with the parameters inside the lasing corridor.
Inside the lasing corridor, a coherent self-oscillation of phonons in the active medium may occur simultaneously with the self-oscillation of the electromagnetic field in the cavity. In Fig. 5a, the region in which coherent phonons are excited is bounded by the solid black curve. In this figure, the stable solutions with coherent phonons belong to the area colored in orange. Usually, the system dynamics is uniquely determined by the external parameters ω a , g, D 0 and is independent on initial conditions of Eqs. (3)-(6), {a, σ , D, b} t=0 . The most intriguing of our results is that there is an area in which two solutions corresponding to different physical regimes are stable. In Fig. 5a, this area lies between the upper dotted line and the dashed line. Moreover, in the area shown in green, the trivial solution and the periodic harmonic solution are stable, while in the area colored in brown, the trivial and periodic non-harmonic solutions are stable. Here, for different values {a, σ , D, b} t=0 , the system arrives at the different stationary state. For example, in the green area, for the initial conditions close to the trivial solution, the system does not begin lasing with time. If the initial conditions are far from the trivial solution, at large times, the conventional laser self-oscillation mode is observed.
It should be noted that each point in Fig. 4 corresponds to a particular system having a specific resonator and active material. The resonator frequency ω a and the Fröhlich interaction constant g are plotted along the axes; the Rabi frequency and the pump rate are fixed. Equations (3)-(6) are solved for each such a system.
The transition from self-oscillations with coherent phonons to chaos. As noted above, in the absence of the electron-phonon interaction, the laser system could be described by the "gas laser" approximation, γ a < γ σ , in which the chaotic regime (the Lorenz attractor) does not arise 29 . Thus, there is no reason to expect a transition of the considered system to the strange attractor regime. www.nature.com/scientificreports/ Computer simulation shows that the electron-phonon interaction significantly changes the laser dynamics. Namely, if the rate of the exciton dephasing γ σ is reduced (note that the inequality γ a < γ σ still holds), then in the region where coherent phonons may arise, areas of chaos appear. Such a decrease in γ σ can be achieved, for example, by lowering the temperature. For such a case, various operational regimes are shown in Fig. 5b.
In Fig. 6, time dynamics |σ (t)| and its spectrum F |σ | (ω) are shown. The parameters g , D 0 , and ω a correspond to the white point in Fig. 5b. With the accuracy of our computer simulation, which is about 10 −7 eV , all frequencies are present in the spectrum. Such a spectrum drastically differs from the discrete spectrum observed for self-oscillations (Fig. 4a). This is an indication of a chaotic regime.

Conclusions
We have considered a conventional single-mode laser consisting of a resonator and a gain medium. Previous studies of conventional laser systems ignored high Raman cross-section of active molecules leading to a strong coupling between excitons and optical phonons. In this work, as far as we know, we consider the effect of the vibrational degrees of freedom on the laser dynamics for the first time.
We describe the exciton-phonon coupling by the Fröhlich Hamiltonian and find that there exists a critical value of the Fröhlich coupling interaction constant g cr . Below g cr , the conventional laser dynamics is observed, above g cr , new operational regimes arise. The value of this constant g cr is about 0.01 eV . This value is realized in Figure 5. Operational regimes of lasing in the parameter plane ω a , g ; (a) γ σ = 0.005eV and (b) γ σ = 0.002 eV . In the areas colored in blue, lasing is absent; the conventional lasing regime occurs in the yellow area; the operational regime with coherent phonons occupies the orange area, the regime of the chaotic operation occurs in the area colored in red; the area where both trivial and common lasing regimes are stable is colored in green (yellow + blue), the area in which both the common laser regime and the regime with coherent phonons are stable is colored in brown (orange + blue); finally, the area in which both trivial and chaotic solutions are stable is colored in purple (red + blue). The other system parameters are D 0 = 0.6 , R = 0.01 eV , γ b = 0.0002 eV , γ a = 0.001 eV , ω b = 0.01 eV , ω σ = 2.4 eV , γ D = 0.001 eV , and γ p = 0.004 eV. www.nature.com/scientificreports/ several types of quantum dots, such as CdSe, CdS, and PbS. The Raman cross-section for these quantum dots is ∼ 4 × 10 −6 Å; it corresponds to the Fröhlich interaction constant g ≃ 0.025 eV , which is of the order of g cr . The parameters ω a , g , and D 0 uniquely determine the type of solution to equations of motion (3)- (6). We have done a comprehensive study of the types of solutions for each set of values of the parameters. Different colored regions in Figs. 1 and 5 correspond to the different types of the obtained solutions.
In one of these regimes, a new type of self-oscillations appears. In this regime, both coherent phonons and the coherent optical field are involved. The period of these self-oscillations is equal to 2π/ω b,gen , where ω b,gen is the frequency of generated coherent phonons,. The spectrum of these oscillations contains a series of equidistant frequencies that differ by the frequency ω b,gen . The transition to the new self-oscillation regime occurs as the Hopf bifurcation, as the pump rate increases. The new regime is characterized by its own pump rate threshold that differs from the lasing threshold.
Another unusual regime is a chaotic laser dynamics, the spectrum of which contains incommensurable frequencies. This regime is realized for the same values of the Fröhlich interaction constant as the previous regime, but for lesser values of the exciton dephasing rate, γ σ ∼ 10 −3 eV , that can be achieved, for example, by lowering temperature.
The most unusual property of the new regimes is that above the second threshold, in certain ranges of the resonator frequency, the pump rate, and the Fröhlich constant, more than one solution of the system equations of motion are stable. In this region, the the initial values of the resonator field, the active medium polarization, the population inversion, and the phonon amplitude govern which regime is realized.
In essence, the new self-oscillations described in the paper represent a generator of coherent optical phonons, which may serve as a basis for phonon nanolaser. Such a device can be used for a substantial resolution enhancement in image processing. In addition, the generation of coherent phonons is important for optomechanical applications 27 , in quantum information technology 44 , and can be used for quantum operations and protocols 45,46 .

Methods
Master equation for the system density matrix. To describe the phenomena under consideration, we use the technique of the Lindblad equation. This equation describes the time-evolution of the density matrix 23,24,37 where the superoperators and describe energy relaxations of the cavity mode and the exciton, respectively 24,37 . The superoperator describes incoherent pumping of the QD exciton 47 . Note that by incoherent pumping we mean the pumping of the high-lying energy levels by a coherent or incoherent light whose frequency is higher than ω σ . Subsequent nonradiative relaxation of electronic states to the exciton state |e� produces the pumping of the exciton.
The superoperator describes the relaxation of the phase of the exciton dipole moment. This superoperator is responsible for dephasing-the phase destruction of the exciton dipole moment. The quantities γ a , γ deph , γ D , and γ p have the meaning of the rates of respective relaxation and pump processes. The last term in Eq. (26) has the form It describes the energy relaxation of optical phonons 24,37 .
dσ/dt = (−iω σ − γ σ /2)σ + i� R aD, www.nature.com/scientificreports/ where the self-oscillation frequency ω gen is determined by the "frequency pulling" formula 23,24,40 : In the general case, this self-oscillating solution becomes unstable if two conditions, γ a > γ D + γ σ and D 0 > D th γ a (γ a + γ D + γ σ )/γ σ (γ a − γ D − γ σ ) , are met. In this case, a regime of deterministic chaos arises, which in the phase space of variables, {a, σ , D} takes the form of the Lorentz attractor 23,49 . But, since we consider the case of a "gas laser" ( γ a ≪ γ σ ) 29 , even the first of the conditions listed above is never satisfied, and there is no reason to expect a transition to the strange attractor regime.
Effective change in the exciton frequency due to a constant displacement of the QD nuclei under incoherent pumping. We begin with considering the QD dynamics in the absence of an external field, including the resonator field. This helps us to interpret the regimes of the QD generation in the cavity more clearly.
In the absence of an external field and a resonator, we can put R = 0 . Equations (3)-(6) are reduced to The trivial stationary solution to this system has the form: Let us consider the dynamics of small deviations of the QD dipole moment σ from this equilibrium position. Substituting σ = σ st + δσ , b = b st + δb , and D = D st + δD into Eqs. (45)-(48), assuming ω σ >> g (that is always satisfied, since ω σ lies in the optical range and g is in the terahertz range), and implying that the value of δb is small in comparison with b st , we obtain where we have introduced the notation for the phonon effective frequency shift of the QD 50 : Thus, the oscillations of the exciton dipole moment near the equilibrium position occur not at the exciton eigenfrequency ω σ but at the shifted frequency ω σ + ω sh . The displacement ω sh is always negative and proportional to the value D 0 + 1 , which has the meaning of the population of the excited state of the exciton at a nonzero pump rate, D 0 > −1. a(t) = γ p + γ D D 0 − D 0,th /2γ a e −iω gen t , σ (t) = i ω gen − ω a − γ a /2 � R γ p + γ D 2γ a D 0 − D 0,th e −iω gen t , ω gen = ω σ γ σ + ω a γ a γ σ + γ a .